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METHOD FOR COMPREHENSIVE IDENTIFICATION OF CELL 
LINEAGE SPECIFIC GENES 
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This application claims priority to U.S. Provisional application no. 
60/440,510, filed on January 16, 2003, the disclosure of which is incorporated herein 
by reference. 

10 FIELD OF THE INVENTION 

The present invention relates generally to the field of identification of 
expressed genes and more particularly to method for characterization of cell lineage 
specific genes. 

15 DISCUSSION OF RELATED ART 

Embryonic and somatic stem cells are considered to offer potential therapy 
for a large spectrum of diseases. Parkinson's, cardiomyopathies, and diabetes are a 
small subset of the potential diseases that could benefit firom the effective 
application of these cells. However, one of the major impediments to the effective 

20 use of embryonic stem cells to treat a broad range of diseases is the lack of sufficient 
knowledge of the mechanisms leading to the differentiation of the required cell 
lineages. 

A cell's lineage defines its relationship to the multipotent precursor that gave 
rise to it. The molecular mechanisms controlling the decision to differentiate 

25 towards one of two or more alternative lineages are of great interest for 

understanding the basic biology of embryonic development as well as homeostasis 
within somatic tissues. Understanding these mechanisms is also important for the 
practical application of stem cell therapy. 

The use of cell-type specific gene expression has historically been a key 

30 molecular method for defining cell types and cell lineage relationships. For 

example, comparison of genes expressed in embryonic stem (ES) cells with those 
expressed in neural stem cells (NSCs) and hematopoietic stem cells (HSCs) suggests 



a closer relationship between ES and NSCs than HSCs since there are more genes 
expressed in common between ES cells and NSCs. More recent application of 
global approaches to surveying gene expression now make it possible to move 
beyond cell type identification to molecular phenotyping based on the expressed 
5 complement of genes. This phenotype can imply function and, to a first 

approximation, cell type is defined at the molecular level by the genes expressed. 

The ability of global approaches for surveying gene expression to give 
insights into function is illustrated by a recent study employing microarray 
technology to provide a comprehensive look at those genes that are expressed within 

10 mouse embryonic stem cells and compare these to genes expressed in neural and 

hematopoietic stem cells ( Ramalho-Santos et al., 2002, Ivanova et al., 2002). The 
results from these studies provide a profile of genes common to the three stem cell 
compartments and these serve as clues to the functions that may be required to 
maintain cells in an imdifferentiated state. Expression profiles from these cells were 

15 compared to those of the lateral ventricle of the brain and the main population of the 
bone marrow and genes showing elevated expression in one or more stem cell 
populations were identified. These studies revealed 216 or 283 genes (Ramalho- 
Santos, 2002 and Ivanova et aL, 2002, respectively) that are elevated in their 
expression in all three stem cell compartments and imply a core set of signaling 

20 pathways that may contribute to maintaining the properties of stem cells. 

Additionally, the observation that there are fewer differences between ES cells and 
NSCs may imply that these cell types are more closely related to each other than to 
HSCs. 

Similar molecular profiles for other multipotent somatic stem cell 
25 compartments and differentiated cell types could reveal much about the functions 
and relationships between these different cell types. Large-scale efforts 
characterizing genes expressed within tissues, often as a function of an additional 
parameter (e.g. age, disease state), have been made by a number of laboratories. A 
public repository of many of these results is available through the NCBI's Gene 
30 Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/entrez/querv.fcgi7db 

=geo) . GEO includes data from large-scale characterizations using Serial Analysis 
of Gene Expression (SAGE) and microarray detection platforms. While these data 



are useful for identifying specific genes within a tissue that responds to an 
experimental variable, the complex cell mixtures within most tissues limit the value 
of such studies for providing a profile of the genes expressed within any specific cell 
type. 

5 Genome wide phenotyping of gene expression within specific cell types can 

be accomplished if a means of isolating the relevant cell type is available. In the 
example cited above, effective methods exist for either growing the cells (ES cells, 
NSCs), or in the case of HSCs isolating relatively pure populations (from an 
unlimited starting cell population, bone marrow, using FAGS). There are, however, 

10 a number of caveats and limitations to applying similar microarray based approaches 
for defining the genes expressed during the differentiation of these cells to specific 
cell lineages. For in vitro differentiation model systems, two means to address this 
issue are: 1) to optimize the firaction of cells that differentiate towards the target 
lineage and/or 2) to physically isolate or select for the desired lineage. In practice, 

15 the first of these methods is generally employed through the inclusion of growth 
factors or ligands, often empirically identified, which can shift the ratio of cells 
differentiating towards a given lineage. While helpful, pure populations of cells 
directed towards only a single target lineage are rarely obtained in practice. For 
many cell lineages that are induced based on cell-cell interactions during the 

20 differentiation process, directing differentiation towards a single lineage may not be 
possible without a-priori knowledge of all of the relevant signals, if at all. 

Physical isolation of target lineages is feasible and presents an efficient 
method for their recovery at higji purity. In the case of the hematopoietic system, 
cell svirface markers allowing identification of most cell compartments are available 

25 and such approaches are already possible in principle. A potentially powerful 
alternative for other cell lineages, where cell surface markers are generally not 
known or available, is the use of FAGS sorted cells to recover cells that are marked 
by EGFP expression within a specific cell lineage. This approach has been applied 
to define the molecular phenotype of Drosophila germline cells by linking EGFP 

30 expression to the vasa gene (Sano et al., 2002), to define oligodendrocytic lineages 
derived from mouse ES cells using an Olig 2-GFP knock-in (Xian et al., 2003). 

While the above approaches are useful for comparative analysis of expressed 
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genes at selected times, many genes, often those of greatest interest, are expressed 
only transiently during the establishment of a cell lineage and will not be reflected in 
the expression profile from difiTerentiated cells. To address this issue, in one 
approach, the fate of cells that express a known marker gene transiently has been 
5 developed (Zinyk et aL, 1998). A bi-genic system in which Cre recombinase 

expression is linked to expression from a gene specific to the precursor cell type is 
present in combination with a reporter that is activated permanently by 
recombination only when Cre is expressed. Rearrangements leading to reporter 
expression then occxu- in the precursor cell and mark its progeny regardless of 

10 whether they continue to express the marker gene that caused the initial expression 
of Cre. This methodology is an effective means of following the fate of a precursor 
cell in the forward direction. In this method, the stem cell carries both a Cre 
recombinase vector expressed under the control of the promoter for a known gene 
and an EGFP (or other) reporter that is activated on Cre expression through a 

1 5: permanent rearrangement within the vector. Hence, transient activation of Cre 

recombinase within a multipotent precursor induces a rearrangement that results in 
EGFP expression from the constitutive Pgk promoter in this cell as well as its 
progeny. Although it is then possible to trace the fate of a precursor cell, there are 
limitations to the utility of this method for determining the molecular profile of 

20 genes expressed within specific cell lineages. Isolation of EGFP expressing cells 

from a differentiating population will allow recovery of the cell and its progeny but 
it will not be possible to assign individual genes to specific cell types within the 
lineage. A second disadvantage is that the precursor specific gene must be known in 
advance. 

25 Currently there are no methodologies available for rapid characterization of 

cell lineage genes to obtain a molecular profile of cells of all or most of the changes 
in gene expression during differentiation. Accordingly, there is an ongoing need to 
develop new approaches to studying cell lineage specific gene expression. 



30 SUMMARY OF THE INVENTION 

The present invention provides methods and compositions for rapid and 
comprehensive identification of cell lineage specific genes. The method of the 
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present invention comprises identification of cells destined toward a particular 
lineage. The steps of the method of the present invention are as follows. The steps 
required for retrospective gene expression analysis are to: 1) establish an embryonic 
stem cell line with a recombinase excision-dependent, cell type specific fluorescent 
protein reporter, 2) use the cell line for recombinase-gene trapping in a library 
fashion where individual cells will not be derived, but rather -10,000 - 30,000 
different insertions will be kept together in a mixed cell culture, 3) isolate cells 
differentiating towards the marked cell lineage on the basis of fluorescent protein 
expression using FACS, 4) recover short sequence tags from the trapped genes in a 
high-throughput fashion and identify the tagged genes. 

The steps of the present invention are accomplished by using a vector system 
comprising the elements for cell lineage targeting, gene trapping and high- 
throughput analysis of the trapped genes. In one embodiment two vectors are used. 
One vector is termed herein as the cell lineage targeting vector comprising a cell 
lineage specific gene, the promoter for that gene, a deletion sequence generally 
targeted by a recombinase, a selectable marker and a reporter gene. The selectable 
marker enables the isolation of cells destined toward the selected lineage. A second 
vector is then introduced into the cells. This vector is termed herein as the gene-trap 
vector. This vector comprises the elements for the modified serial analysis of gene 
expression (MAGE), a recombinase and a selectable marker. 

In another embodiment, this invention provides a vector system comprising 
the cell lineage targeting vector and the gene-trap vector for identification and 
temporal characterization of cell lineage specific genes. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a schematic diagram of retrospective gene expression analysis for 
identifying genes expressed within the lineage of a differentiated cell. 

Figure 2 is a schematic diagram of the expected gene identifications based on 
profiling differentiated cells, cells marked by the currently used forward approach 
and the reverse gene expression analysis of the present invention. A is a 
Constitutive gene, B is specific to I, C is specific to II, D is specific to la, E is 
specific to lb, F is specific to Ha, G is specific to lib, H is specific to Ibl and I is 



specific to Ib2. 

Figure 3 is a schematic representation of the cell lineage targeting vector. 

Figure 4 is a schematic representation of the gene-trap vector. 

Figure 5 is a representation of the modified SAGE approach for primary 
5 reporter identification from gene trap cell lines. 

Figure 6 is a representation of the gene-trap vector structure and integration 
of this vector into an endogenous gene. 

Figure 7 is a representation of the structure of a Cre gene-trap vecotr 
following integration into an endogenous gene. 
10 Figure 8 is a representation of the structxire of a targeting vector for 

integration of a Cre-dependent Emerald reporter construct into the 3' non-translated 
region of the Synapsis I gene. 

Figure 9 is a representation of fluorescence showing lineage marking with a 
Cre-gene trap vector with or without Tomaxifen. 
15 Figure lOa-h is a representation of the fluorescence of humanized EGFP and 

Emerald fluorescent proteins when expressed fir-om the CMV promoter in 293 cells. 
Figure 10a and 10b are maps of plasmids encoding the humanized GFP (pTRGS- 
green plasmid) and Emerald (pcDNA3 -Emerald plasmid) fluorescent proteins. A 
histogram the fluorescence intensity is shown in Figures 10c and lOd; scatter plots 
20 are shown in Figures lOe and 1 Of and the fluorescence and phase contrast images are 
shown in Figures lOg and lOh. 

Figure 1 1 is a representation of illustrative results from the application of 
MAGE to a pool ofgene trap marked cell lines. Electrophoresed PGR products are 
shown in Panel A. Electrophoresis of all PGR products in Panel A except the 232 
25 bp is shown in Panel B. 

Figure 12 is a representation of an illustrative plasmid usefiil for inverse 

PGR. 

DETAILED DESCRIPTION OF THE INVENTION 
30 The method of the present invention is based on the following concept. It is 

desirable to associate a specific differentiated cell type with a molecular profile of 
all of the changes in gene expression that occiuxed during its derivation. A 

-6- 



conceptual scheme allowing such a profile to be obtained is shown in Figure 1 . In 
this figure, a simple lineage relationship is shown in which there is a stem cell that 
can give rise to two multipotent precursor cells I and II. These in turn can lead to 
the differentiated cell types la and lb or Ila and lib. Genes are represented by the 
5 letters A - G, where all cells express A, B is specific to multipotent precursor I, C is 
specific to multipotent precursor II and D - G are lineage specific markers for the 
differentiated cells la - lib as shown in the figure, a molecular profile of the 
differentiated cell la would be constituted by the expression of genes A, B, and D; 
reflecting the derivation of la firom I. 

10 In the method of the present invention, a reporter gene expression is driven 

under control of the differentiated cell specific gene D promoter. Hence, the only 
cells that will express the reporter gene are the differentiated cells la and if the 
reporter gene is a fluorescent protein, these can be isolated by FACS. However, la 
cells will only express EGFP if a rearrangement in the reporter construct has 

15 occurred due to the expression of a recombinase. This occurs only if the gene 

promoter driving the recombinase was active at some point in the derivation of cell 
la. Further, if recombinase is randomly integrated into different genes in a pool of 
the starting stem cell population (i.e. in a library fashion), different genes can be 
sampled for their expression during the derivation of cell type la. 

20 A comparison of the information generated by expression profiling methods 

currently available and the method of the present invention is shown in Figure 2. 
The advantage of the retrospective gene expression approach is that it makes it 
possible to start with a well defined differentiated cell type in which genes have 
been trapped randomly with a gene-trap vector and in which a cell lineage gene is 

25 marked, and look backwards into the gene expression history of just the lineage 
giving rise to that specific cell type. By this method a through analysis of genes 
expressed during differentiation all the way back to the point where the genes were 
trapped, can be done Further, it will capture genes expressed within the lineage even 
if their expression is only transient and they are no longer expressed in the 

30 differentiated cell. If retrospective gene expression data are generated for two or 
more lineages derived from a common stem cell, it will also be possible to 
reconstruct the molecular relationship between the lineages from gene expression 



data obtained by the method of the present invention. 

Accordingly, the present invention provides a method for rapid identification 
of cell lineage specific genes in embryonic stem cells. The method of the present 
invention comprises introducing into a stem cell line a cell lineage specific gene 
5 whose expression, detected by a reporter protein, is dependent upon a recombinase 
excision event. The recombinase is randomly integrated into different genes. Also 
present in the vector carrying the recombinase gene are elements of a high 
throughput analysis for identification of trapped genes. When cells destined toward 
the selected cell lineage are identified (based on the expression of the reporter 
10 protein), the trapped genes (indicating genes expressed at some point during 
differentiation) can be identified in a high-throughput fashion. 

No previous methodology can comprehensively define the history of gene 
expression, exclusive of those expressed constitutively, within a specific cell 
lineage. Applied on a broader scale, knowing the history of gene expression for cell 

15 types that are derived from a common progenitor will allow establishment of the 
point at which they diverged and the molecular phenotype of the precursor cell 
since, at some point in their history, a common set of genes will be identified for the 
different cell types. This makes it possible to know the molecular details of the 
relationships between all of the cell types formed following differentiation of ES 

20 cells in vitro. 

The ability to define genes expressed throughout the history of a given cell 
lineage would have enormous potential utility. There are nxraierous applications for 
this information that would facilitate the objective of defining cell lineages for 
transplantation. For example, from this set of genes, it is possible to cull those that 

25 would provide useful markers for various stages of differentiation towards a desired 
lineage. The ability to mark precursors at various stages in their differentiation 
towards different cell lineages for recovery, e.g. with EGFP, makes it possible to test 
their ability to contribute to a desired tissue. Some of the genes are likely to be 
potential regulatory molecules. It is possible to use these to direct differentiation to 

30 a desired lineage. Cell surface molecules that are specific to early stages of the 
desired lineage may also be identified. This would offer the possibility of raising 



-8- 



antibodies to be used in recovering such cells in the absence of genetic marking and 
may be particularly important for application to transplantation of such cells in 
humans. Identification of receptors or components of signal transduction pathways 
could offer new insights into biologically active agents that may facilitate 
5 differentiation of cells towards specific cell lineages. The same set of molecules 
may be pharmacological targets for modulating the fimction of stem or progenitor 
cells. 

The present method can also be used in the context of a whole animal. 
Animals carrying the necessary recombinase dependent reporter protein marked cell 

10 type specific gene can be generated (e.g. firom the aMHC knock in) as follows. 
Bone marrow recovered firom these cells can be treated essentially as described 
herein for the stem cells and the cells can then be retumed to the animal. 
Altematively retrovimses are highly effective in transducing cells in vivo. Hence, at . 
least in those cases where sufficient numbers of stem cells are present for targeting 

15 and sufficient nimibers of the specific cell lineage to be assessed can be recovered, 
the gene expression history of the cell type from its stem cell progenitor can be 
traced in'an adult animal. For example, for use in vivo, the Cre-HSVtk recombinase 
may be advantageous since it would allow elimination of constitutively expressed 
trapped genes in situ, prior to the start of the gene to lineage linking protocol, 

20 through the administration of gancyclovir. 

Further, this method can be used for investigating the divergence of genes by 
identification of temporally expressed genes in two or more lineages. In siunmary, 
this technology can be used for fully defining the relationships between cell lineage 
and gene expression in vitro and, potentially, in vivo. 

25 The elements required for the present invention i.e., cell lineage targeting, 

gene-trapping and high throughput analysis can be introduced in one or more vectors 
into cells. In one embodiment, the elements are introduced as two separate vectors. 
The vectors can be introduced in the cell together or in any sequential order. One 
vector is termed as the cell lineage targeting vector. This vector comprises a cell 

30 lineage gene promoter (non-targeted vector for non-targeted recombination), a 
selectable marker and a reporter gene. Further, a pair of sequences targeted by a 
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recombinase are present flanking a polyadenylation site. In one embodiment, the 
cell lineage targeting vector comprises a cell lineage specific gene, a selectable 
marker, a reporter gene and the recombinase target sites flanking the 
polyadenylation site (targeted vector for targeted recombination). In one 
5 embodiment, the selectable marker may be present in the sequence flanked by the 
recombinase target sites. An example of a targeted vector is shown in Figure 3. A 
second vector is termed as the gene-trap vector. It comprises the gene for a 
recombinase (such as Cre or Flp), a selectable marker and elements for carrying out 
a rapid identification of trapped genes by the modified serial analysis of gene 
10 expression (MAGE). An example of such a vector is shown in Figure 4. 

Any cell lineage specific gene may be used to constmct the cell specific 
lineage targeting vector. Such cell specific genes are well known in the art. For 
example, to identify the lineage of cardiomyocytes, alpha MHC gene and its 
promoter can be used. For identifying the lineage of neurons, neuron specific genes 
15 such as synapsis or neuron specific enoleise can be used. For identifying glial cells, 
glial fibrillary acidic protein (GFAP) gene can be used. Other examples of cell 
specific lineiage genes are available on the NCBI-GEO web cite which is easily 
accessible and well known to those skilled in the art. 

The selectable marker in the cell specific targeting vector can be based on 
20 any positive or negative selection approach. The selectable marker facilitates the 

selection of host cells transformed or transfected with the vector. Positive selection 
marker genes are those that encode a protein such as G418, hygromycin, puromycin 
and blastocidin, which confer resistance to certain drugs and proteins allowing for 
positive selection. The negative selection includes conditional negative selection, 
25 such as HSV-tk in the presence of gancyclovir or acyclovir, and nitroreductase in the 
presence of metronidezole 

The reporter gene in the cell specific targeting vector is any sequence of 
DNA that encodes for a protein which is detectable by an assay. In a preferred 
embodiment, the protein is a fluorescent protein. For example, the green 
30 fluorescence protein gene (GFP) isolated fi-om the jellyfish Aequorea Victoria, has 
become available as a reporter in prokaryotes and eukaryotes. The gfp gene encodes 
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a protein which fluoresces when excited by violet or blue-green light. Variants of 
GFP are also available. One such variant is the enhanced GFP or EGFP (which is 
shown in Figure 1 as being regulated by the promoter of the cell specific lineage 
gene, oMHC. Another example of a reporter protein is red fluorescent protein 
5 (RFP) and yellow fluorescent protein (YFP). Further, proteins which are detectable 
via histochemical stains (e.g., b-galactosidase, alkaline phosphatase) can be used. 
Expression of the reporter gene is dependent upon the promoter of the cell lineage 
specific gene and on the deletion of a recombinase target sequence. In the absence 
of a recombinase, the recombinase target sequence is not deleted and the 

10 transcription of the reporter gene is prevented by the polyA site. The recombinase 
target sequence depends upon the recombinase in the gene-trap vector. Thus, if the 
recombinase is Cre, the recombinase target sequence comprises the loxP sites and if 
the recombinase is FLP, the recombinase target sequence comprises the fil 
sequences. Further, fiisions between two protein that confer the functions of each 

1 5 may also be used (such as b-GEO). It will be recognized by those skilled in the art 
that reporter proteins other than fluorescent proteins can be alternatively used. Thus, 
any reporter protein which allows the isolation of cells in which rearrangement by 
gene-trapping has taken place, can be used. Such reporter proteins include magnetic 
tags, positive selection such as puromycin, lacZ protein and the like. 

20 In the gene trap vector is present the DNA sequence encoding a recombinase. 

Examples of recombinase known in the art include Cre and FLP. In one 
embodiment, the recombinase may be fiised to fluorescent protein such as HCRed, 
or Thymidine Kinase or nitroreductase. In the gene trap vector are also present 
elements for MAGE (described in U.S. patent application serial no. 10/227,719, 

25 filed on 8/26/02, incorporated herein by reference) and as described herein. 

The vectors of the present invention can be introduced into suitable host cells 
by standard methods of transfection including lipofection, precipitation, infection, 
electroporation, microinjection and the like. Such methods are well known in the art 
(see for example, see Sambrook et al.. Molecular Cloning A Laboratory Manual, 
30 2nd edition. Cold Spring Harbor Laboratory Press, 2001 

In one illustration, the steps of the present method are as follows. A cell 
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lineage specific gene (such as oMHC as an indicator for cardiomyocytes or synapsin 
I as an indicator of neuronal cells) that can be used for marking lineage is selected. 
A recombinase (such as Cre or flp) dependent reporter gene is then knocked in to 
mark the expression of the selected cell lineage gene. The sequence deleted by the 
5 recombinase also comprises a selectable marker such as a PGK driven neomycin 

resistance gene. An example of such a vector is shown in figure 3. By the use of this 
vector, cells destined toward the selected lineage are identified by the selectable 
marker and these cells can be isolated. 

In the next step, cells are tranfected with a gene-trap vector comprising a 
10 recombinase. The gene-trap vector also comprises elements for a high throughput 
detection assay. Preferably cells without base level recombinase are selected. In 
one embodiment, an inducible recombinase (such as tamoxifen dependent Cre 
recombinase) can be used. An example of a gene-trap vector is shown in Figure 4. 

The expressed genes at various stages of differentiation can be characterized 
15 by using the modified serial analysis of gene expression (MAGE) technique. Key to 
the ability of the present technology to provide a global profile of lineage specific 
gene expression is a means of identifying gene trap insertion sites efficiently. 
Elements allowing high efficiency acquisition of sequence tags that identify such 
sites are incorporated within the splice jxmctions of the Cre-Gene Trap Vector to be 
20 used. It is this feature of the gene trap vector that permits a modified version of the 
Serial Amplification of Gene Expression (SAGE, Velculescu et al., 1995) 
technology to be utilized in identification of trapped genes in a high throughput 
format. This technology is referred to as MAGE and is described in detail below. 

The Modified SAGE technology (MAGE) is a high throughput method for 
25 the identification of sequence tags resulting firom gene trap vector integration events. 
The basis of this technology is shown in Figure 5. The first element on which it 
depends is the incorporation of type IIS endonuclease restriction (such as Bsgl and 
Bpml) recognition sequences adjacent to the splice acceptor and splice donor 
elements within the HTP gene trap vector. These type IIS restriction endonucleases 
30 have the property that each cleaves the DNA at a position 16 nucleotides adjacent to 
the recognition sequence where the composition of the 16 nucleotides is irrelevant. 
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Other examples of type IIS sites are BsmFI, Mmel and Fokl. Using the example of 
Bsgl and Bpml, as shown in Figure 6, we have taken advantage of this property to 
allow the amplification of either 15 or 14 nucleotides of the endogenous gene 
sequence adjacent to the SA and SD elements of the gene trap vector, respectively. 
5 This allows differential amplification of endogenous gene sequence from cDNAs to 
messages that result from transcripts initiating from the endogenous gene promoter 
when Bsgl is used or the Pgk promoter when Bpml is used. Hence, in the case of 
Bsgl, the resulting products will reflect the relative expression level from the 
marked gene when assaying mixed pools while Bpml will result in relatively even 

10 levels of amplification products. This becomes important below. 

Following this, bits of unknown sequence information can be identified 
because these are separated by repeats of a known sequence. In the present 
application, this is accomplished by ligating the PGR products with the aid of a 
restriction endonuclease cleavage site present in both the imiversal primer and 

15 adjacent vector sequence. The ligated strings of sequence tags are then cloned and 
sequenced. Thus, each sequence tag representing each member present in a pool of 
marked genes regardless of the absolute expression level, can be sequenced. Since 
transcripts expressed from the Pgk promoter will be present at relatively equal 
levels, use of the SD junction fragments is optimal. Data similar to the following 

20 sequence string can be obtained: 

5TCTAGACAGTCTGGA GAGNNNNNNNNNN]^^ 

GAGNNNNNNNNNNNNN^ 

NNTCTAGA NNNNTSTNNNNN^ 

GAGNNNNN3* (SEQ ID NO: 1). Each repeating unit is 32 nucleotides long and 
25 contains 16 nucleotides that are derived from a discrete gene trap event (the splice 

donor AG plus 14 as underlined) and can be used to identify the insertion site. 

Inversion of the repeats is possible; however, this event is easily recognized. 

In one embodiment, to address the issue of different RNAs being expressed 

at different levels which can skew representation of certain RNAs, instead of RT 
30 PGR, inverse PGR can be used. Inverse polymerase chain reaction (IPGR) is an 

extension of the polymerase chain reaction that permits the amplification of regions 

that flank any DNA segment of known sequence, either upstream or downstream 
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(see U.S. Pat. No. 4,994,370, specifically incorporated herein by reference in its 
entirety). The essence of IPCR is that, by circularizing a restriction enzyme fi-agment 
containing a region of known sequence plus flanking DNA, PGR can be performed 
using oligonucleotides whose sequence is taken from the single region of known 
5 sequence and oriented with respect to one another such that their 5* to 3' extension 
products proceed toward each other by going "around the circle" through what 
originally was flanking DNA. This leads to the amplification of DNA strands 
containing what was originally flanking DNA. The advantage of a technique such as 
IPCR, with respect to the current invention, is that using a single primer set one may 
10 amplify a representative sample of insertion jimctions from a particular group of 
individuals. 

An illustrative example of inverse PGR is provided below. An illustrative 
plasmid is shown in Figure 12. For inverse PGR, genomic DNA is pooled fi-om the 
cells expressing the reporter gene. An aliquot is digested with Mspl The digested 
1 5 DNA is ligatged to circularize and amplified using an LTR primer and biotinylated 
MmeI-1 primer. The amplicon is immobilized on a strepavidin tube and digested 
with Mmel to expose genomic tags. The tags are ligated to a universal oligo, 
amplified and purified by HPLG. The tags can then be concatamerized and ligated 
into plasmid vectors. 

20 In one embodiment, the cloning vector containing two different cloning ends 

(in one embodiment. Sad and NotI) are used. The tag concatamers are derived by 
ligation of the tags at a dilute concentration in the presence of 0.5:1 Molar ratio of 
each of two non-phosphorylated DNA adapter molecules (in this embodiment the 
adapters would contain both Xbal overhangs for ligation to the tag cloning ends, 

25 while on the opposite end one adapter would have a NotI overhang and the other 
would have a Sad overhang). The use of non-phosrphorylated adapters minimizes 
the formation of any inhibitory side reaction products. The use of different ends on 
the cloning vector allows the more efficient directional cloning of the concatamers. 
The adapters are added to the reaction before addition of the tag monomers in order 

30 to maximize monomer tag addition while simultaneously minimizing self ligation of 
monomers into mini DNA circles (which would result in a loss of critical material). 
The ligation of an adapter to one end of a concatamer prevents circularization of that 
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molecule while allowing continued addition of monomer tags to the other end until 
that too becomes ligated to an adapter. Those concatamers which have different 
overhangs at each end (e. g. Sad at one end and NotI at the other) will clone with 
high efficiency into the recipient vector. 
5 This invention is further described in the Examples presented below which 

are intended to be illustrative and not restrictive. 

EXAMPLE 1 

In a further illustration of this invention, the method can be carried out in 
10 three steps. Initially, an embryonic stem (ES) cell is created in which a specific cell 
lineage is marked by knocking in a Cre-mediated recombination dei)endent EGFP 
vector. Second, a highthroughput (HTP)-gene-trap vector that both marks genes 
with a tamoxifen dependent Cre recombinase and allows detection of expression of 
the marked gene by expression of hcRFP (or altematively allows negative selection 
15 of the gene by expression of herpes simplex vims thymidine kinase), is created. 
Finally, the ES cell line is utilized in conjimction with the gene trap vector as 
follows to identify genes expressed at different stages in the differentiation towards 
the marked cell lineage. 

Following transfection of the ES cell line with the gene trap vector and 
20 selection for its integration, cells can be first sorted for lack of RFP expression, 
induced to differentiate, and treated with tamoxifen at various times. Due to Cre 
mediated re-arrangement of the EGFP reporter, cells that carry a gene trap event in a 
gene that is expressed in the marked cell lineage will, ultimately, express EGFP. 
These cells can be recovered by FACS, and the genes incorporating the gene trap 
25 vector can be determined using a HTP-sequence tag identification technique 
(MAGE). 

1) ES cell line derivation. To establish the initial knock in ES cell line a 
vector of the stmcture shown in Figure 3 and allowing Cre dependent activation of 
EGFP expression firom a cell lineage specific gene is created. The example shown 
30 in Figure 3 uses the aMHC promoter which is useful for marking cardiomyocytes. 



- 15- 



The construct can be electroporated into W4 ES cells and colonies can be selected in 
G418 and scored for correct integration at both the 5' and 3' ends. Sufficient 
colonies are scored to establish between 3 and S independent lines. One or more 
lines are selected for further use on the basis of efficient differentiation towards 
5 cardiomyocyte lineages. 

2) Gene trap vector construction. A gene trap vector incorporating the 
sequences necessary for both HTP sequence tag acquisition by MAGE and reverse 
orientation retroviral packaging can be constructed as shown in Figure 4. This 
vector carries a tamoxifen dependent Cre recombinase fused to the fluorescent 

10 protein HCRedl. Hence, on integration into an endogenous gene, expression fi-om 
the gene will both lead to Cre expression and be detectable through fluorescence in 
the Red chaimel. The vector is packaged and utilized to transduce the aMHC 
conditional EGFP marked cell line derived as described above. Cells are selected in 
the presence of puromycin and the resulting colonies are pooled. Approximately 

15 10,000 to 30,000 colonies are desired, which using an appropriately tittered viral 

stock, are achievable using approximately 10x15 cm culture dishes. Colonies can 
be pooled and used in the protocol described below. 

3) Identification of lineage specific genes. The third and informative 
component of this technology is to recover lineage specific (in the current example 

20 cardiomyocyte) genes fi-om the pool of gene trapped cells. This is accomplished as 
follows. First, cells are sorted for those that do not express detectable levels of RFP. 
It is important here that cells that are negative for RFP do not express Cre 
recombinase at functional levels. Hence, a sample of the non-RFP expressing 
population is grown in the presence of tamoxifen (for 48 hours) and assessed for 

25 recombination at the aMHC conditional EGFP gene by PCR. In the event that 

significant levels of recombination are occurring, it can be concluded that RFP is not 
a sufficiently sensitive marker and the gene trap vector to incorporate and Cre- 
HSVtk fusion can be reconstructed. In this case cells can be selected against HSVtk 
expression using gancyclovir and the efficacy of this selection can be confirmed by 

30 analysis of recombination at the aMHC conditional EGFP gene. If HSVtk is not a 
sufficiently sensitive negative selection, progressively disabled versions of Cre 
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recombinase through sub-optimal codon usage substitutions can be created as is 
known to those skilled in the art. The effect of the negative selection at this step is 
to remove from the gene trap marked cell population those genes that are active in 
iminduced cells. 

5 A second level of characterization of the non-expressing cells is performed to 

define the fiill range of tagged genes. For this analysis MAGE is performed from 
the splice donor side of the vector on the pool of un-induced non-RFP expressing 
cells. If we have tagged the desired 30,000 cells, there should be a corresponding 
30,000 sequence tags to acquire or 30,000 x 32 bp of non-redxmdant sequence. This 

10 is 960,000 bp. Each sequence run is expected to yield about 500 bp of information 
so approximately 1,900 sequence reactions or 20 microtitre plates of sequence will 
be required. This can be performed on the MegaBase capillary sequencer within 3 - 
4 days. For approximately 3 X coverage to detect 95% of the tagged genes, the 
estimated time to complete the sequencing is about 3 weeks. Establishing the fiill 

15 complement of genes that could be detected provides an important base line against 
which to measure the subset that is detected in the desired cell lineage as described 
below. 

Following sorting for cells that do not express RFP (or, altematively, grow in 
the presence of gancyclovir), these cells will be induced towards cardiomyocyte 

20 Uneages using an optimized induction regime. Tamoxifen is introduced into the 
culture media at progressively later stages of induction in different cultures where, 
initially, selected intervals (such as 1 day intervals) can be used. In each culture, 
cells are treated with tamoxifen (e.g., for 24 hours). During this treatment, Cre 
expression from active genes will allow recombination at the cell lineage specific 

25 conditional EGFP gene. Recombination will occur in all cells that carry a marked 
gene that is active during the time tamoxifen is added regardless of whether they are 
destine to become, in the example here, caridomyocytes or alternative lineages or 
whether the gene remains active even within the cardiomyocyte lineage. Regardless 
of the time of tamoxifen addition, cells will be cultured through 8 days at which 

30 point activation of the aMHC gene should occur within all cells in the 

cardiomyocyte lineage. Only those cells that expressed Cre from a gene that was 
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marked by a gene trap event and expressed in the cardiomyocyte lineage during the 
time tamoxifen was active will activate EGFP expression. Cells are then trypsinized 
and sorted using FACS to recover these cells. Sequence tags from the marked genes 
will be identified using MAGE from the splice donor Bpml site. These sequence 
tags will identify sets of genes expressed specifically in the cardiomyocyte lineage at 
various points in its derivation from ES cells. If sufficiently large numbers of gene 
trap events are scored, this method can be used to define and temporally order the 
expression of the large majority of genes that are specific for the cardiomyocyte 
lineage. 

EXAMPLE 2 

This example describes the construction of another cell lineage specific 
targeting vector. In this example, a Cre-recombinase excision dependent Synapsin I- 
specific Emerald reporter construct is described. The neural specific gene, 
Synapsin I, is known to be expressed in neurons derived from retinoic acid (RA) 
treated ES cell cultures (Finley et al; 1996). The targeting construct shown in Figure 
7 was prepared and an Asel to Pvul fragment was isolated and electroporated into 
the W4 ES cell line (Taconics). 

G418 resistant colonies can be amplified and assessed for correct targeting 
into the Synapsin I gene by BamHl and Xbal digestion and Southem blotting to 
assess correct targeting at the 5* and 3* ends respectively. In this illustration, the 
targeting construct comprised the fluorescent protein Emerald. The expression of 
this Emerald from the Synapsin I promoter following Cre recombinase mediated 
excision of the triple polyA stop signal (or lox-STOP-lox cassette, abbreviated 
XTX) allows the identification and isolation of these cells. 

EXAMPLES 

This example describes the Construction and transduction efficiency of a 
lenti-viral based gene trap vector. An example is shown in Figure 6. Elements 
required for gene-trapping fimctions include a reporter gene (here EGFP) 
downstream of a splice acceptor such that, on integration into an intron of an 
endogenous gene, the reporter will become spliced into the endogenous message 
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allowing its expression. In most cases this also disrupts function of the endogenous 
gene. An internal ribosome entry site (IRES) is placed 5' to the EGFP sequence to 
allow its expression regardless of the reading frame of the endogenous transcript. 
To insure that ribosomes initiating from the endogenous transcript start codon will 
5 not occlude the IRES or result in a fusion protein, a series of translation termination 
codons are placed 5' to the IRES. The vector also carries a selectable marker (here 
neo) driven from a constitutive promoter (Pgk) and followed by a splice donor to 
allow selection of stably transfected cell lines on integration into an endogenous 
gene. Not shown are viral packaging sites that allow reverse orientation packaging 
10 into a self-inactivating (SIN) lenti virus. 

Elements allowing high efficiency acquisition of sequence tags are 
incorporated within the splice junctions. This is a key feature of the vector that 
permits this technology to be utilized in identification of trapped genes in a high 
throughput format. 

15 In this embodiment, the following modifications were made (shown in figure 

8). First, the fluorescent protein reporter has been substituted for an - 
HSVtk/CreERT2 cassette, which encodes a fusion protein between the herpes 
simplex virus thymidine kinase gene and a tamoxifen dependent Cre recombinase 
(Feil et aL, 1997); Although only the Cre recombinase component is required for 

20 the core gene-trap methodology, variations of the method may take advantage of 
either the HSVtk or tamoxifen dependence of Cre in this fusion protein and these 
variations are discussed in the experimental methods section. Second, to allow use 
in cells that already carry a neo resistance gene, the neo gene has been replaced with 
a puromycin resistance gene (PURO) in the present construct. Other elements, 

25 including the SA and SD sequences modified for higji-throughput sequence tag 
analysis are maintained in the Cre gene trap vector. 

The function of both the Cre-recombinase and the loxStoplox (XTX) 
elements used in the Cre-gene-trap-vector and the Cre-dependent-Synapsin-I- 
Emerald-reporter constmcts, respectively, were tested as shown in Figure 9. Here, 

30 the Pgk promoter was used to drive expression from either the gene trap vector as 
shown in Figure 8 or an EGFP reporter containing the XTX element used in the 
Synapsin construct as shown in Figure 7. NIH 293 cells were tranfected with either 
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the lineage marking vector alone or this vector plus the Cre gene trap vector and 
grown in either the presence or absence of tamoxifen (Figure 8). No expression of 
EGFP from the lineage-marking vector was foimd in the absence of the Cre gene- 
trap vector. In the presence of the Cre gene trap vector and tamoxifen, but not in the 
5 absence of tamoxifen, robust EGFP expression was observed in co-transfected cells. 

EXAMPLE 4 

This example describes the isolation of fluorescent protein expressing cells 
by FACS. An important component of this invention is the ability to recover 

10 reporter gene expressing cells in the absence of significant levels of non-expressing 
cells. As an illustration of this embodiment, a comparison of the fluorescence of 
humanized GFP and Emerald fluorescent proteins when expressed from the CMV 
promoter in HEK293 cells was carried out. The plasmids used are shown in Figure 
10a and 10b. The results are shown in Figure lOc-h. Figures 10c,e and g are hGFP 

15 while Figures lOd, f and h are Emerald. HEK293 cells were transiently transfected 
with each plasmid using lipofectamine 2000. 36 hours later cells were trypinized 
and sorted on a BD FACS Vantage instrument. Figures 10c and d are histograms of 
fluorescence intensity (FLl) on a log scale; Figures lOe and f are Scatter plot log 
fluorescence against forward scatter (FSC);. Figures lOg and h are fluorescence 

20 images of transfectants obtained just prior to sorting using a Nikon Eclipse TE300 
fluorescent microscope and a SPOT diagnostics camera. Images were processed 
using SPOT diagnostics software version 3.0.4. 

EXAMPLES 

25 This example describes a high-throughput vector dependent sequence tag 

recovery. As proof of principle, the generation and identification of Sequence Tags 
from trapped genes by MAGE from the splice donor jimctions present in a small 
pool of PI 9 EC cell gene trap lines was performed. RNAs were isolated using 
GIT/phenol extraction and polyadenylated messages were selected on oligo dT 

30 cellulose by standard methods. First strand cDNA synthesis primed with oligo dT 
was performed using superscript II (Invitrogen) using standard conditions. A 
control sample in which reverse transcriptase was omitted was also prepared. RNA 
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was hydrolyzed using NaOH. NaOH was neutralized and cDNAs were recovered by 
ethanol precipitation. Second strand synthesis was rimed using Biotinylated 
neotop2 primer (5'-B-CCGCTTTTCTGGATTCAT-3' - SEQ ED NO:2) and 
extended using the large fragment of E. coli DNA polymerase. Double stranded 
5 cDNA was degiested with Bpml and incubated with streptavidin coated PGR tubes 
for 3 minutes at 37C. Following binding tubes were washed times with 150ul of 
150mM NaCl in TE (10 mM Tris-HCl, pH 7.5, 1 mM EDTA). The MAGE 
universal adapter (5'-TCTAGAGGACTGCGTGGGCGA-3' - (SEQ LD NO:3); 5'- 
CCTCGCCCACGCAGTCCTCTAGANN-3' -(SEQ ID NO:4) (16.6 nmoles) was 

10 added in 50 ul of ligation buffer plus 2 ul of T4 ligase and tubes and tubes were 

incubated for 2 hours at 15C. Following 2 washes as previously, 1 ul of 50 uM of 
each MAGE PGR primer (5'-CCTCGCCCACGCAGTCCTC-3' (SEQ ID NO:5), 
5'CGGCTGGGTGTGGCGGAC-3' - (SEQ ID NO:6)) was added in 100 ul of 
Platinum Taq (Invirogen) PGR reaction buffer containing 0.2 mM of each of dATP, 

15 dGTP, dCTP and dTTP, 2 mMMgG12 and 0.5 units of Platinum Taq polymerase. 
Thermal cycling was performed where 35 cycles of 94C for 0.75 minutes, 60C for 
0.75 minutes and 72C for 0.75 minutes. Samples of the resulting PGR products 
were electrophoresed on an agarose gel as shown in Panel A of Figure 1 1. The 
predicted band of 232 bp is present in the lanes in which reverse transcriptase was 

20 present in the initial cDNA synthesis reaction but is absent from lanes whre reverse 
transcriptase was initted. The remaining PGR products were pooled, digested with 
Xbal and electrophoresed on an 8% polyacrylamide as shown in Panel B. The 
predicted 32 nt fragment containing the sequence tags was isolated and incubated in 
20 ul of T4 DNA Ugation buffer for between 20 and 60 minutes and used for 

25 transformation of competent HB 101. Transformed colonies of HB 101 were selected 
and used in preparation of DNA for sequencing by standard techniques. The data 
from the sequencing is shown in Table 1. 



Table 1 



Sequence Tag 


Gene 


GGGGAGAGGGGGAGGT 


erythroid differentiation regulator 
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(SEQ ID NO:7) 




AvjvjvjAvjvjAvjA IKjx J\Kj 1 

(SEQ ID NO:8) 




AU 1 AUA XK^L 1^1^ 1 1 

(SEQ ID NO:9) 


unigcne ciusier jj^inv^i 


OCjCj 1 ILr i Ci 1 U AU 1 CI 

(SEQIDNOrlO) 


unigene cluster, LAr-rciaicu rcirovirai 
elements 


CCC AOOA i Cj 1 A 1 ACjC 1 
(SEQ ID NO: 11) 


IMAvjii : 34 1 o3 y o 


AOl C 1 CjVjAAACC lOl C 
(SEQIDNO:12) 


cJJJNA clone xlJ 12Jr 1 1 J 


AC 1 ACAl C 1 CC 1 ACC 1 
(SEQ ID NO: 13) 


CLilNA Clone VJ3U143i5UO J 


ALrC i AALf 1 CrC 1 LrC 1 LrC 

(SEQIDNO:14) 


mouse i^oi, ui-JVi-r>riz.J aog-c-UD-u- 
Ul.sl 


ACjCj 1 AAo 1 AC AAvjC 1 0 
(SEQ ID NO: 15) 


mouse Hoi ujoUDiu.xi 


ACCTCACjTTCjA 1 1 CC 1 
(SEQ ID NO: 16) 


CUNA clone A43UU30rU3 J 


AOCTTTCjAAT 1 C A 1 LrA 
(SEQ ID NO: 17) 


, nspo4-l 


1 1 ICICAOCjCtI AUCCl 
(SEQIDNO:18) 


albumin 


1 CCLjC 1 CAA 1 Cj 1 ACC 1 

(SEQIDNO:19) 


unknown 


ACCC 1 1 CjCtAC 1 CAACtCf 
(SEQ ID NO:20) 


unknown 


(SEQIDNO:21) 


UIUvIllJ WXl 


cacgggggcggagcct 

(SEQ ID NO:22) 


unknown 
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(SEQ ID NO:23) 


unKnown 


(SEQ ID NO:24) 


iitilmnwn 

UXIAXIV/ WIX 


CGCCATCGCTCTAGCT 
(SEQ ED NO:25) 


unknown 



Sequencing revealed that the concatamers consisted of the predicted 
vector/universal primer sequences separated by 16 nucleotide long tags. Blast 
searches of the tags revealed seven unknown sequences (i.e. not present in the NCBI 
mouse EST or non-redimdant sequence databases) and twelve known sequences 
comprising predicted exons from albumin, HSP84, actin binding protein, erythroid 
differentiation regulatory protein and others. 

EXAMPLE 6 

In another illustration of the embodiment, further sequence tag sequences 
were generated as described in Example 5. The sequences are provided below in 
Table 2. 



Table 2 



well 
no. 


Sequence Tag 


Gene 


A2 


AGGCTCATCAGCTGAC 
(SEQIDNO:26) 


RIKEN cDNA 1600014E20 (chr. 2, 2F2) 
embryo E6.5-E8.5, placenta 


A4 


AGGTGTGGATCCAGAG 
(SEQIDNO:27) 


RIKEN cDNA 9130023F12 gene (chr. 
11) mammary gland; Tcell; head; urinary 
bladder; thymus; lung; brain; spontaneous 
tumor, metastatic to mammary. Stem cell 
origin; neural retina 


A7 


AGGTAGGTCTCGATCT 
(SEQ ID NO:28) 


no identification 
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A Q 


Ar:A ATA r*r; AT A PPPAri 

(SEQ ID NO:29) 


no lucniiiicaiion 


D / 


A nnnr* a TTTr^T A tt a 0 
(SEQ ID NO:30) 


no lacnuncaiion 




A r^r^T A A TC^ A f^r^ a a a i^r^ 

ALri^ 1 AA i AL/V^ AAAVJ\^ 

(SEQIDNO:31) 


3G2) numerous tissues 


C6 


AGGATCCACTGTTTAC 
(SEQ ID NO:32) 


no identification 


D12 


AGGTCTGGAGTAACCA 
CACATAGATGTTAGTT 

A A/^A/^A/^A A A A i^T^ A A 

AACj ALf AO AAAACj i AA 
CTGGAGACTTCCTCAC 

A AT^T-An 
AAl vjAvj 

(SEQ ID NO:33) 


partial Bpml digest product, 
phospholipase c-like 2 gene, opp ori 
relative to gene, prob cryptic splice (chr. 
9) 


ElO 


AGGTCCTGCCTCAGCA 


repetitive 


F2 


AGGACTGATTGTGGTG 


G protein-coupled receptor 39 (Gpr39, 
cnr. 1, iCiO) nean, lesiis, emoryo, leius, 
whole brain; visual cortex* 




ACjAIAACjI 1 ILrl IdCj 
(SEQIDNO:36) 


no identification 


"DO 


ALrVj 1 ALr 1 AvjVjVjAI^U 
(SEQIDNO:37) 


no laenuncanon 


00 


a r;r;T A TO A Trt A p A or^T 
(SEQ ID NO:38) 


no lu.cnuiiL'duun 


H3 


AGGTGTACGAATGCGA 
(SEQ ID NO:39) 


no identification 



*Lineage restricted gene trapped by XTX - Gpr39 

From the foregoing, it will be obvious to those skilled in the art that various 
modifications in the methods described herein can be made without departing from 
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the spirit and scope of the invention. Accordingly, the invention may be embodied 
in other specific forms without departing from the essential characteristics thereof. 
The embodiments and examples presented herein are therefore to be considered as 
illustrative and not restrictive. 
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